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Abstract 

The field equations of a generalized f(R) type gravity model, in which there is an arbitrary 
coupling between matter and geometry, are obtained. The equations of motion for test particles 
are derived from a variational principle in the particular case in which the Lagrange density of the 
matter is an arbitrary function of the energy-density of the matter only. Generally, the motion 
is non- geodesic, and takes place in the presence of an extra force orthogonal to the four- velocity. 
The Newtonian limit of the model is also considered. The perihelion precession of an elliptical 
planetary orbit in the presence of an extra force is obtained in a general form, and the magnitude 
of the extra gravitational effects is constrained in the case of a constant extra force by using Solar 
System observations. 
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I. INTRODUCTION 



Recent astrophysical observations have provided the astonishing result that around 95- 
96% of the content of the Universe is in the form of dark matter + dark energy, with only 
about 4-5% being represented by baryonic matter [l]. More intriguingly, around 70% of 
the energy-density is in the form of what is called " dark energy" , and is responsible for the 
acceleration of the distant type la supernovae [2J. Hence, today's models of astrophysics 
and cosmology face two fundamental problems, the dark energy problem, and the dark 
matter problem, respectively. Although in recent years many different suggestions have 
been proposed to overcome these issues, a satisfactory answer has yet to be obtained. 

A very promising way to explain the observational data is to assume that at large scales 
the Einstein gravity model of general relativity breaks down, and a more general action 
describes the gravitational field. Theoretical models in which the standard Einstein-Hilbert 
action is replaced by an arbitrary function of the Ricci scalar R, first proposed in [3], have 
been extensively investigated lately. Cosmic acceleration can be explained by f(R) gravity 
4], and the conditions of viable cosmological models have been derived in jr|. In the context 
of the Solar System regime, severe weak field constraints seem to rule out most of the models 
proposed so ta B, B, althou^ viable mo de,s do exist fi fl U ~ 



11] . The possibility that 
'or dark 



the galactic dynamic of massive test particles can be understood without the neec 
matter was also considered in the framework of f(R) gravity models 



For a review of f{R) generalized gravity models see 17]. 



131, 



14, 



15 



16] . 



A generalization of the f{R) gravity theories was proposed in 18] by including in the 
theory an explicit coupling of an arbitrary function of the Ricci scalar R with the matter 
Lagrangian density L m . As a result of the coupling the motion of the massive particles is 
non-geodesic, and an extra force, orthogonal to the four-velocity, arises. The connections 
with MOND and the Pioneer anomaly were also explored. The implications of the non- 
minimal coupling on the stellar equilibrium were investigated in [jj]], where constraints on 
the coupling were also obtained. An inequality which expresses a necessary and sufficient 
condition to avoid the Dolgov-Kawasaki instability for the model was derived in 2(|. The 
relation between the model with geometry-matter coupling and ordinary scalar-tensor grav- 
ity, or scalar-tensor theories which include non-standard couplings between the scalar and 
matter was studied in 21[ . The motion of matter in such theories as well as the dark matter 



2 



problem in galaxies was also analyzed. In the specific case where both the action and the 
coupling are linear in R the action leads to a theory of gravity which includes higher order 
derivatives of the matter fields without introducing more dynamics in the gravity sector 



221 ] . The equivalence between a scalar theory and the model with the non-minimal coupling 



of the scalar curvature and matter was considered in [23j. This equivalence allows for the 
calculation of the PPN parameters (3 and 7, which may lead to a better understanding of the 
weak-field limit of f(R) theories. Different forms for the matter Lagrangian density £ m , and 
the resulting extra-force, were considered in 24(, and it was shown that more natural forms 
for C m do not imply the vanishing of the extra-force. The impact on the classical equiva- 
lence between different Lagrangian descriptions of a perfect fluid was also analyzed. Similar 
couplings between gravitation and matter have also been considered as possible explanations 



25]. 



for the accelerated expansion of the universe and of the dark energy in 

In all the previous studies of the models with matter-geometry coupling, the matter part 
in the coupling was represented by the Lagrangian density of the matter, while the geometric 
part was considered to be an arbitrary function of the Ricci scalar. We may call this class 
of models as modified gravity models with linear matter-geometry coupling. However, more 
general models, in which the matter part in the coupling is an arbitrary function of the 
Lagrangian density of the matter, can also be constructed, and they represent the natural 
generalization of the models with linear matter coupling. It is the purpose of this Letter 
to present the field equations of a generalized gravity model, in which there is an arbitrary 
coupling between matter and geometry. In this class of models the energy-momentum tensor 
of the matter is generally not conserved. The equations of motion of the test particles are 
also obtained, by using a variational principle, in the particular case in which the Lagrange 
function of the matter is a function of the density of the matter only. In this model the 
motion is non-geodesic. The study of the Newtonian limit shows that the matter-geometry 
coupling induces a supplementary acceleration of the test particles. This acceleration may 
be responsible for the constancy of the galactic rotation curves, which is usually attributed 
to the presence of the dark matter. As an observational test of the model we consider the 
perihelion precession of an elliptic planetary orbit in the presence of the extra force. The 
precession angle is derived in a general form, and from the observed value of the perihelion 
precession of the planet Mercury the magnitude of the extra acceleration is constrained. 
The present Letter is organized as follows. The field equations of the model are derived 
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in Section II. The equations of motion of the test particles and their Newtonian limit are 
considered in Section III for a particular model in which the matter Lagrangian is a function 
of the density only. We discuss and conclude our results in Section IV. 



II. GRAVITATIONAL FIELD EQUATIONS IN f(R) TYPE MODELS WITH AR- 
BITRARY COUPLING BETWEEN MATTER AND GEOMETRY 

The most general action for a f(R) type modified gravity involving an arbitrary coupling 
between matter and geometry is given, in a system of units with 8nG = c = 1, by 

=gd A x, (1) 



S 



l -h{R) + G{L m )f 2 {R) 



where fi(R), % = 1,2 are arbitrary functions of the Ricci scalar R, while G (L m ) is an 
arbitrary function of the matter Lagrangian density L m . The only requirement for the 
functions fi, i — 1,2 and G is to be analytical function of the Ricci scalar R and L m , 
respectively, that is, they must possess a Taylor series expansion about any point. When 
fi(R) = R, f2(R) = 1 an d G (L m ) = L m , we recover standard general relativity. With 
f2(R) = 1 and G (L m ) = L m we obtain the f(R) generalized gravity models. The case 
G (L m ) = 1 + XL m , X = constant, corresponds to the (linear) coupling between matter and 
geometry, considered in [3] - j^]. 

We define the energy-momentum tensor of the matter as 

By assuming that the Lagrangian density L m of the matter depends only on the metric 
tensor components, and not on its derivatives, we obtain T^ u = L m g^ u — 2dL m /dg^ u . 

Varying the action with respect to the metric tensor g^ v we obtain the field equations of 
the model as 



F 1 (R)R^ U - -h(R)g^ + (g^D - V^V,) F 1 (R) = -2G (L m ) F 2 {R)R^ U 



1 

T 

-2(g^D-V ll V u )G(L m )F 2 (R)- 



f 2 (R) [K (L m ) L m -G (L m )) 9flv + f 2 (R)K (L m ) T^, (3) 

where we denoted Fi(R) = df i {R)/dR ) i = 1,2 and K (L m ) = dG (L m ) /dL m , respectively. 
For G (L m ) = L m and by rescaling the function f 2 (R) so that f 2 (R) — > 1 + Xf 2 (R), we 
reobtain the field equations proposed in 



By contracting the field equations given by Eq. we obtain the scalar equation 

3D [Fx (A) + 2G (L m ) F 2 (R)] + [F^R) + 2G (L m ) F 2 (R)] R - 

2h{R) + Af 2 (R) [K (L m ) L m -G (L m )] = K (L m ) f 2 (R)T, (4) 

where T = T£. By taking the covariant divergence of Eq.Q, with the use of the mathe- 
matical identity V [a'(R)R llv - a(R)g^/2 + {g^U - V M V^) a(R)] = [26], where a(R) is 
an arbitrary function of the Ricci scalar R and a'(R) = da/dR, we obtain 

dL 

= In [f 2 (R)K {L m )\ {L m9fiu - T^} = 2V^ In [fo{R)K {L m )\ j^. (5) 

The requirement of the conservation of the energy-momentum tensor of matter, V^T^j, = 
0, gives an effective functional relation between the matter Lagrangian density and the 
functions /^(-R) and K (L m ), 

BT 

Vln[/ 2 (^(L m )]^=0. (6) 

Thus, once the matter Lagrangian density is known, by an appropriate choice of the 
functions G (L m ) and f 2 {R) one can construct, at least in principle, conservative models 
with arbitrary matter-geometry coupling. 



III. MODELS WITH ARBITRARY DENSITY-DEPENDENT MATTER LA- 
GRANGIAN 

As a specific case of generalized gravity models with arbitrary matter-geometry coupling, 
we consider the case in which the matter Lagrangian density is an arbitrary function of the 
energy density of the matter p only, so that L m = L m (p) . Then the energy-momentum 
tensor of the matter is given by 

where the four-velocity w M = dx^/ds satisfies the condition g^u^u^ = 1, and we have also 
used the relation 5p = (1/2) p (g^ v — u^Uy) 5g^ u . 

The energy-momentum tensor given by Eq. ([7j) can be written in a form similar to the 
perfect fluid case if we assume that the thermodynamic pressure p obeys a barotropic equa- 
tion of state, so that p — p(p). The perfect-fluid type representation can be obtained by 
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assuming that the matter Lagrangian satisfies the equations 

P^^ = P + pTl(p) + p(p), (8) 

and 

P^^-L m (p)=p(p) } (9) 

respectively, where II (p) is an arbitrary function of the density Substituting the term 
pdL m (p) I dp from Eq. (jSJ) into Eq. IQ gives the matter Lagrangian as L m (p) — p + pU (p). 
With this form of L m , Eq. ([8]) gives the following differential equation for II (p), 

P—^=P(J>)> ( 10 ) 



with the general solution given by 

Jo P Jo P P 

Therefore the matter Lagrangian and the energy-momentum tensor can be written as 

L m (p)=p(l + fj^j -p(p), (12) 

and 

IT = [p + p(p)+ P U (p)] uV - p (p) gr, (13) 

respectively. From a physical point of view II (p) can be interpreted as the elastic (defor- 
mation) potential energy of the body, and therefore Eq. (fT3"j) corresponds to the energy- 
momentum tensor of a compressible elastic isotropic system. 

By imposing the condition of the conservation of the matter current, V u {pu v ) = 0, and 
with the use of the identity u v V v u^ = d 2 x^ jds 2 + T^ x u u u x , from Eq. (jSJ) we obtain the 
equation of motion of a test particle in the modified gravity model as 

^+i>v = r, (M) 

where 

/" = -V, In j/ a (J2)ff [L m (p)] (uV - <T) ■ (15) 

The extra-force / M , generated due to the presence of the coupling between matter and 
geometry, is perpendicular to the four-velocity, f^u^ = 0. The equation of motion Eq. ffT4"|) 
can be obtained from the variational principle 

5S P = S J L p ds = 5 J y/Q y /g /lv uf 1 u u d,s = 0, (16) 



where S p and L p = \fQ^fg^uFvF are the action and the Lagrangian density for the test 
particles, respectively, and 

Q = h(R)K[L m (p)}^^. (17) 



To prove this result we start with the Lagrange equations corresponding to the action ( 1161) , 

d fdL p \ dL p _ o (ig) 



ds \du x J dx x 

Since dL p /du x = y/Qu x and dL p /dx x = (1/2) y/Qg^xvPu" + (1/2) Q,\/Q, a straightfor- 
ward calculation gives the equations of motion of the particle as 

d 2 x^ 

~~dl? 



+ T^u u u x + [u»u v - g^) V u In VQ = 0. (19) 



By simple identification with the equation of motion of the modified gravity model with 
arbitrary matter-geometry coupling, given by Eq. (jl4j) . we obtain the explicit form of \/Q 
as given by Eq. (TlT|) . 

The variational principle ffl6l) can be used to study the Newtonian limit of the model. 



In the limit of weak gravitational fields, ds ~ yjl + 20 — iPdt ~ (1 + (f) — v 1 /2) dt, where 
is the Newtonian potential and v is the usual tridimensional velocity of the particle. By 
representing the function y/Q as 

v/Q = f 2 (R)K [L m (p)} = 1 + U (r } L m (p) , dJ ^) , (20) 

where U << 1, the equations of motion of the particle can be obtained from the variational 
principle 



S 

and are given by 



dL m (p)\ v 



U[R,L m (p),^^j +0- - 



dt = 0, (21) 



-V0 -VU = a N + a E , (22) 



where 3n = — V0 is the usual Newtonian gravitational acceleration, and ag = — V£7 is a 
supplementary acceleration induced due to the coupling between matter and geometry. 

An estimation of the effect of the extra-force, generated by the coupling between matter 
and geometry, on the orbital parameters of the motion of the planets around the Sun can 
be obtained in a simple way by using the properties of the Runge-Lenz vector, defined as 
A — v xL — ae r , where v is the velocity relative to the Sun, with mass M , of a planet 



of mass m, f = re r is the two-body position vector, p = /j,v is the relative momentum, 
p, = mM Q / (m + -Me) is the reduced mass, L = r xp = [ir 2 9k is the angular momentum, 



and a = GmM P) 



271 ] . For an elliptical orbit of eccentricity e, major semi-axis a, and period 



T, the equation of the orbit is given by (L 2 / fia) r^ 1 = 1 + ecos#. The Runge-Lenz vector 
can be expressed as A = (L 2 / fir — a J e r — rLee, and its derivative with respect to the polar 
angle 9 is given by dA/d9 = r 2 [dV(r)/dr — a/r 2 ]ee, where V(r) is the potential of the 



central force 



271 ] . The potential term consists of the Post- Newtonian potential, Vpjv( 



—a/r — 3a 2 /mr 2 , plus the contribution from the general coupling between matter and 
geometry. Thus we have dA/d9 = r 2 [6a 2 /mr 3 + ma E (r)} eg, where we have also assumed 
that // ~ m. The change in direction A0 of the perihelion with a change of 9 of 2tt is 



obtained as A</> = (1/ae) J Q 27r L x dA/d9 



d9, and it is given by 



AJ nA -U a \ 2 1 L (l-e 2 f 2 [ a E \L 2 (I + e cos 6) /ma] n in /ooN 

A0 = 24tt 3 - - + — - Q V- / — " ^r 1 L cos6d6, 23 

v \TJ 1-e 2 8^me ( a /T) 3 J (1 + ecos^) 2 1 ' 

where we have used the relation a/L = 2n (a/T) /yl — e 2 . The first term of this equation 
corresponds to the standard general relativistic precession of the perihelion of the planets, 
while the second term gives the contribution to the perihelion precession due to the presence 
of the coupling between matter and geometry. 

As an example of the application of Eq. ff23]) we consider the case for which the extra- 
force may be considered as a constant, a E ~ constant, an approximation that could be valid 
for small regions of the space-time. With the use of Eq. ( 1231) one finds for the perihelion 
precession the expression 



6nGM Q 2vraVl - e 2 
a(l -e 2 ) + GM & 



^ = : n 1 + —7^7r *e> (24) 



where we have also used Kepler's third law, T 2 = 4n 2 a 3 /GM & . For the planet Mercury 
a = 57.91 x 10 11 cm, and e = 0.205615, respectively, while M = 1.989 x 10 33 g. With 
these numerical values the first term in Eq. (I1H]) gives the standard general relativistic value 
for the precession angle, (A(p) GR = 42.962 arcsec per century, while the observed value of 
the precession is (A(f)) obs = 43.11 ± 0.21 arcsec per century [28(. Therefore the difference 
(Acj)) E = (A(f>) obs — (A(j)) GR = 0.17 arcsec per century can be attributed to other physical 
effects. Hence the observational constraints requires that the value of the constant a E must 
satisfy the condition a E < 1.28 x 10~ 9 cm/s 2 . This value of a E , obtained from the solar 
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system observations, is somewhat smaller than the value of the extra-acceleration ao ~ 10 8 
cm/s 2 , necessary to explain the "dark matter" properties, as well as the Pioneer anomaly 



181 ] . However, it does not rule out the possibility of the presence of some extra gravitational 
effects acting at both the solar system and galactic levels, since the assumption of a constant 
extra-force is not correct on larger astronomical scales. 

IV. CONCLUSIONS 

In the present Letter we have considered a generalized gravity model with an arbitrary 
coupling between matter and geometry, described by the product of an arbitrary function 
of the Lagrange density of the matter, and an arbitrary function of the Ricci scalar. The 
proposed action represents the most general extension of the standard Hilbert action for the 
gravitational field, S = J [R/2 + L m ] y/ ' —gd A x. The equations of motion corresponding to 
this model show the presence of an extra-force acting on test particles, and the motion is 
generally non-geodesic. The physical implications of such a force have been already analyzed 
in the framework of the generalized gravity model with linear coupling between matter and 



geometry, considered in 18j , and the possible implications for the dark matter problem and 
for the explanation of the Pioneer anomaly have also been investigated. On the other hand, 
the field equations Eqs. (j3J) are equivalent to the Einstein equations of the f(R) model in 
empty space-time, but differ from them, as well as from standard general relativity, in the 
presence of matter. Therefore the predictions of the present model could lead to some major 
differences, as compared to the predictions of standard general relativity, in several problems 
of current interest, like cosmology, gravitational collapse or the generation of gravitational 
waves. The study of these phenomena may also provide some specific signatures and effects, 
which could distinguish and discriminate between the various gravity models. 
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